Neutron Dripline in Odd and Even Mass Calcium and Nickel Nuclei 
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Neutron rich Ca and Ni nuclei have been studied in spherical Relativistic Mean Field formalism 
in co-ordinate space. A delta interaction has been has been adopted to treat the pairing correlations 
for the neutrons. Odd nuclei have been treated in the blocking approximation. The effect of the 
positive energy continuum and the role of pairing in the stability of nuclei have been investigated 
using the resonant-BCS (rBCS) approach. In Ca isotopes, N = 50 is no longer a magic number 
while in Ni nuclei, a new magic number emerges at N = 70. There is a remarkable difference in the 
relative positions of the drip lines for odd and even isotopes. In Ca isotopes, the last bound even and 
odd nuclei are found to be 72 Ca and 59 Ca, respectively. In Ni isotopes, the corresponding nuclei are 
98 Ni and 97 Ni, respectively. The origin of this difference in relative positions of the dripline in even 
and odd isotopes in the two chain is traced to the difference in the single particle level structures 
and consequent modification in the magic numbers in the two elements. Pairing interaction is seen 
to play a major role. The effect of the width of the resonance states on pairing has also been 
investigated. 

PACS numbers: 21.60. Jz,27.40.+z,27.50.+e,27.60.+j 
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■ In recent years, it has been possible to populate and study a number of neutron rich nuclei from fusion evaporation 
reactions using radioactive ion beams as well as from fission fragment studies. The last bound neutrons in such a 
nucleus may lie very close to the continuum and the effect of the positive energy continuum on the structure of such 
nuclei should be studied carefully. Another very important aspect of nuclei near the drip line is the additional stability 
provided by the pairing interaction. It is important to study the effect of pairing by studying the odd mass nuclei 
along with the even mass ones. Although a number of such studies has been undertaken in lighter mass regions, nuclei 
in the medium and the heavy mass regions have not yet been studied in sufficient detail. In the present calculation, 
neutron rich Ca and Ni nuclei have been studied using the Relativistic Mean Field (RMF) formalism in co-ordinate 
ly-j , space. 

RMF theory is a major tool of nuclear structure physics. Very often, the RMF equations are solved by expanding in 
a harmonic oscillator basis. However, it is well known that the basis expansion approach using the harmonic oscillator 
basis cannot explain the density in halo nuclei near the drip line because of slow convergence in the asymptotic region. 
Nonrelativistic Hartree Fock Bogoliubov (HFB) and Relativistic Hartree Bogoliubov (RHB) methods in co-ordinate 
space have emerged as two very accurate approaches of treating the nuclei very close to the drip line. 

The effect of the states in the continuum has been incorporated in most of the calculations by solving the equations 
in coordinate space using the box normalization condition, thus replacing the continuum with a set of discrete positive 
energy states. However, in this case, the single particle energy levels depend on the size of the box chosen, leading to a 
rather unsatisfactory scenario. Nonrelativistic mean field equations involving continuum states have been solved with 
exact boundary conditions QJ3 for zero range and finite range pairing forces. RMF equations involving continuum 
states have also been solvedjjjQ with exact boundary conditions. All these calculations have taken into account 
the effect of the width of the continuum levels also. For example, Cao and Ma Q have also compared their results 
with calculations assuming the continuum levels to be of zero width. We have also used the RMF calculation in 
co-ordinate space including the continuum states to study neutron rich even-even C and Be nuclei in detailQ. It has 
been observed that the results of this method and the more complicated RHB approach are in good agreement 0, H • 
In the present work we study the structure of neutron rich Ca and Ni nuclei beyond 48 Ca and 68 Ni, respectively, 
using the RMF formalism in co-ordinate space with exact boundary conditions, paying particular attention to the 
odd mass nuclei. There have been numerous calculations on even-even neutron rich Ca and Ni nuclei up to the 
neutron drip line using relativistic and nonrelativistic mean field approaches. We list only a few of them which have 
been published recently Nonrelativistic Hartree Fock and Hartree Fock Bogoliubov approach have been utilized by 
Fayans et aZ.0| and Im et ai.0 to study Ca isotopes. Relativistic density dependent Hartree Fock approach Q as 
well as Relativistic Continuum Hartree Bogoliubov method (RCHB)|9j have also been utilized to study neutron rich 
isotopes of Ca and Ni0, |n|. Yadav et al.\lj^ have performed RMF+BCS calculation for Ca and Ni isotopes by 
discretizing the continuum with box normalization. The problem of wrong asymptotic behaviour in the harmonic 
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oscillator basis has been overcome using expansion in a Woods Saxon basis in a RMF calculation [13|. The RCHB 
approach particularly gives a very good description of the nuclei. Terasaki et al.\T^ have studied even-even Ni isotopes 
from the proton to the neutron drip line using HFB approach in three dimensions. However, we have come across no 
mean field calculation for odd Ni isotopes near the drip line. Even in Ca isotopes, the odd mass nuclei have not been 
studied in detail. The aim of the present work is to investigate the odd mass neutron rich Ni and Ca nuclei up to the 
drip line. 

The relativistic mean field theory is well known and readers are referred to Refs. f° r details. We have used the 
force NLSH 16|. We have assumed spherical symmetry for all the nuclei. The BCS calculation has been performed 
adopting a delta interaction to treat the pairing correlations between neutrons, i.e. V = VoS(ri — r 2 ). The usual 
BCS equations now contain contributions from the bound states as well as the resonant continuum. The equations 
involving these states have already been obtained[I|,|3| and have been referred to as resonant-BCS (rBCS) equations. 
We have also included the effect of the width of the positive energy levels. These equations have been solved in the 
co-ordinate basis on a grid of size 0.08 fm. The positive energy resonance solutions are obtained using the scattering 
approach. All the negative energy states beyond N=20 and the positive energy states up to N=82 for which resonance 
solutions are available have been included in the rBCS calculation for studying the Ca isotopes. For the Ni isotopes, 
the corresponding numbers are N=28 and N=126, respectively. We have assumed that beyond 20 fm, the effect of 
nuclear interaction vanishes. We choose Vq= -350 MeV for the strength of the delta-interaction. This value of the 
strength has been used by Yadav et al.\T^ m Ca nuclei. In odd nucleus the last odd nucleon breaks the time reversal 
symmetry of the system. However, it is well known that the bulk quantities like binding energy or radii are not 
affected by the breaking of the symmetry. Since we are interested only in these quantities, we have used the blocking 
approximation to study odd mass nuclei. 

We have calculated the binding energy corresponding to the different levels in odd nuclei. Although the ground 
state spin parity is unknown in most of the nuclei, we find that, except for 75 Ni, the spin parity has been correctly 
predicted wherever the ground state is unambiguously known. In Fig. 1 we plot the binding energy values of the 
all the nuclei studied in the present work which are stable against neutron emission. The experimental or empirical 
values are from Ref. ^tJ- O ne can see that the theoretical values agree with experiment. In Figs. 2 and 3, we have 
plotted the one neutron separation energies for Ca and Ni isotopes, respectively. The one neutron separation energy 
has been defined as S n (Z,N) = B.E.(Z, N) — B.E.(Z,N — 1). In Ni nuclei, the theoretical results are very close to 
experimental measurements. For the Ca isotopes, the agreement is slightly poorer. The calculated results for two 
neutron separation energy for the even mass nuclei are very close to the RCHB results obtained by Zhang et a/.|ll|. 



One of the most interesting differences between the two chains studied is the difference in the positions of the 
dripline for even and odd isotopes. For Ca isotopes, the last odd mass nucleus stable against neutron emission is 59 Ca 
while the corresponding even mass isotope is 72 Ca. On the other hand, the corresponding nuclei in Ni are 97 Ni and 
98 Ni, respectively. In other words, all the Ni isotopes are stable against neutron emission up to N = 70 beyond which 
neither the even nor the odd isotopes are stable for this element. For Ca the odd isotopes are stable below N = 40 
while stable even isotopes range up to N = 52. 

To understand this dramatic difference we next study the single particle neutron levels. In Fig. 4 we have plotted 
the single particle energy levels near the Fermi level for even-even Ca nuclei. The Fermi level is indicated by the solid 
line. Near the drip line the Fermi level rises slowly and finally becomes positive beyond Ca. The states l<?g/2 and 
2^5/2 start as positive energy states but as the neutron number increases they become loosely bound states. Because 
of the absence of the centrifugal barrier, the level 3si/ 2 does not have a resonant solution. The levels 3si/ 2 , 2d 5 / 2 
and 2d3/ 2 lie very close to each other and are close to the l<7g/ 2 level in 66-72 Ca. The shell closure is taken to be 
the neutron number for which the pairing energy vanishes. We find that in Ca isotopes, no shell closure is observed 
beyond N=40 because of the modification in the energy levels. The Fermi level beyond N = 40 increases very slowly 
and lies very close to the continuum. The single particle density near is large near it. Hence pairing correlations which 
bind the loosely bound nucleons can develop in even-even nuclei. As pointed out by Meng et al.\tf$ these conditions 
indicate that a neutron halo should develop in these nuclei. In Fig. 5, we have plotted the neutron radii of the nuclei 
studied. The sudden increase in radius is clearly observed. Importantly, in Ca isotopes, the level I59/2 lies either 
in the continuum or just bound around N=40. Hence the stability of the isotopes depends crucially on the pairing 
interaction and the contribution of the states in the continuum. Thus the odd mass Ca isotopes beyond N = 40 are 
unstable against neutron emission. There is a similar situation in light nuclei. The dripline nuclei n Li and 1 Be are 
known to have two neutron halos. It is known that neither 10 Li nor 13 Be are bound. The pairing interaction stabilizes 
the last neutron pair so that the above two dripline nuclei are stable against neutron emission. 

On the other hand, the Ni isotopes present a completely different picture. The single particle levels are shown in 
Fig. 6. Here, the Fermi level is deeper and abruptly becomes positive beyond A=98. The levels 3si/ 2 , I57/2, 2d 3 / 2 , 
and 2d 5 / 2 he very close to each other. All of them except the first start as positive energy states and become bound 
states at higher neutron number. The intruder orbit l/in/2 lies much higher in energy. This is a consequence of the 
quenching of the spin orbit splitting in neutron rich nuclei which has been obtained in many mean field calculations 
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and also experimentally observed in light nuclei. Thus N = 70 becomes a new magic number. Because the occupied 
levels are much deeper, one does not expect a halo like structure and the binding of the odd nucleus in this case 
does not depend overmuch on the pairing. Hence we obtain both even and odd mass nuclei which are stable against 
particle emission up to the new magic number N=70. The calculated neutron radius values are shown in Fig. 5. One 
can see that the neutron radius shows a shell closure at N — 50, and more importantly, at N = 70. 

It is interesting to study the effect of the width of the resonant states in the nuclei very close to the drip line. In 
Table I, we present the results for some of these nuclei. The width of different levels, the binding energies and the 
neutron radii (r n ) for 70 ' 72 Ca and 96 > 98 Ni are presented. Pairing correlations usually become stronger if the width of 
the resonances be nee lectedHQ thus increasing the total binding energy and decreasing the neutron radius. One can 
see that as 98 Ni is a closed shell nuclei, there is no effect of the resonance states beyond N = 70 in the binding energy 
and the radius. 

To summarize, neutron rich Ca and Ni nuclei have been studied using RMF formalism in the co-ordinate space. 
Exact boundary conditions have been used to obtain the wave function of the states in the continuum. A delta inter- 
action has been used for interaction between neutrons. Odd nuclei have been treated in the blocking approximation. 
The study of odd mass neutron rich nuclei in this mass region reveals interesting features. There is a remarkable 
difference in the relative positions of the drip lines for odd and even isotopes in the two elements. In Ca isotopes, the 
last bound even and odd nuclei are found to be Ca and 59 Ca, respectively. In Ni isotopes, the corresponding nuclei 
are 98 Ni and 97 Ni, respectively. The origin of this difference in the relative positions of the drip line in even and odd 
isotopes in the two chains is traced to the difference in the single particle level structures and the modification in 
the magic numbers in the two elements. In Ca isotopes, N — 50 is no longer a magic number while in Ni nuclei, a 
new magic number emerges at N = 70. The level density near the continuum is higher in Ca isotopes and pairing 
interaction plays a crucial role in the their stability. In odd mass nuclei, the pairing interaction is weaker and not able 
to bind the neutrons in the halo orbitals. If one neglects the width of the resonance states, the pairing correlation 
increases. 
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TABLE I: Width of different resonant levels, and binding energies and neutron radii for calculations involving width (RMFW) 
and ignoring width (RMFN). 



Nucleus 


level 


Width 
MeV 


B.E.(MeV) 
RMFW RMFN 


rn(fm) 
RMFW RMFN 


vu Ca 


2d 5/2 

lg7/2 

lhn/2 
2d 3/2 


0.102 
0.874 
1.524 
0.593 


468.30 468.35 


4.728 4.581 


"Ca 


2d 5/2 

lg7/2 

lhn/2 
2d 3/2 


0.013 
0.448 
1.501 
0.478 


468.35 468.47 


4.994 4.710 


ao Ni 


lhn/2 
lh 9 /2 

1113/2 
2 7 /2 


0.033 
1.678 
2.474 
2.289 


665.02 665.32 


4.952 4.926 




lhn/2 
lh 9 /2 

lil3/2 
2f7/2 


0.018 
1.403 
2.446 
2.000 


667.07 667.07 


4.965 4.965 
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FIG. 1: Experimental and calculated binding energy per nucleon in Ca and Ni nuclei. 




FIG. 2: Neutron separation energy in Ca isotopes. 
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FIG. 3: Neutron separation energy in Ni isotopes. 
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FIG. 4: Single particle neutron energy levels in even Ca isotopes. The Fermi level is shown by the continuous line. 
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FIG. 5: Neutron radius in neutron rich Ca and Ni nuclei. 
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FIG. 6: Single particle neutron energy levels in even Ni isotopes. The Fermi level is shown by the continuous line. 



